Surface criticality at a dynamic phase transition 



(N 

o 



c3 
c3 



o 
o 



>: 

00, 
(TV 

in 

co; 

ON' 

o: 



X 



Hyunhang Park and Michel Pleimling 1 

'Department of Physics, Virginia Tech, Blacksburg, VA 24060-0435, USA 

In order to elucidate the role of surfaces at nonequilibrium phase transitions we consider kinetic 
Ising models with surfaces subjected to a periodic oscillating magnetic field. Whereas the corre- 
sponding bulk system undergoes a continuous nonequilibrium phase transition characterized by the 
exponents of the equilibrium Ising model, we find that the nonequilibrium surface exponents do not 
coincide with those of the equilibrium critical surface. In addition, in three space dimensions the 
surface phase diagram of the nonequilibrium system differs markedly from that of the equilibrium 
system. 
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The ubiquity of nonequilibrium steady states in nature 
constitutes a permanent reminder of the challenges en- 
countered when trying to understand interacting many- 
body systems far from equilibrium. Whereas in some 
instances, as for example paradigmatic transport models 
[l[ or driven diffusive systems notable progress has 
been achieved in understanding nonequilibrium steady 
states, a common theoretical framework remains elusive. 
This is especially true in cases where steady states are in- 
fluenced by the presence of surfaces or interfaces, which 
can change properties even deep inside the bulk [3l4l7|. 

Nonequilibrium phase transitions form an interesting 
class of phenomena that share many commonalities with 
their equilibrium counterparts. For example, for con- 
tinuous transitions different universality classes, char- 
acterized by different sets of critical exponents, have 
been identified. Well-known examples can be found in 
driven diffusive systems Q, at absorbing phase transi- 
tions [l8|, [l9[ or in magnetic systems subjected to a pe- 
riodically oscillating external field [2(1 EHl . For some ab- 
sorbing phase transitions, as for example directed perco- 
lation, the surface critical properties have been studied 
to some extent, see Q and references therein. 

Kinetic ferromagnets in a periodically oscillating mag- 
netic field display as a function of the field frequency a 
nonequilibrium phase transition between a dynamically 
disordered phase at low frequencies and a dynamically 
ordered phase at high frequencies. Let us assume that 
the magnetization is aligned with the direction of the ex- 
ternal field. If the field now reverses direction, the system 
becomes metastable and tries to reverse its magnetization 
through the nucleation of droplets that are aligned with 
the field. If the period of the field is large compared to 
the metastable lifetime, then the metastable state com- 
pletely decays before the field reverses direction again, 
i.e. the ferromagnet is able to 'follow' the field, yielding 
a time-dependent magnetization that oscillates symmet- 
rically about zero. The dynamically ordered phase is 
obtained when the period of the field is small compared 
to the metastable lifetime, thus that the system is not 
able to fully decay from the metastable state before the 
field changes direction again. The magnetization then os- 



cillates about a non-zero value. This behavior has been 
studied theoretically in a large range of systems, as for 
example the Isin g |22h24|| , Heisenberg 25|, l26j , or Blume- 
Emery-Griffiths [27( models, too name but a few. Possi- 
ble experimental realizations have been discussed in Co 
films on Cu(OOl) [28} as well as in [Co/Pt]3 magnetic 
multilayers [29j. Of special interest in the following is 
the kinetic Ising model in an oscillating field that is dis- 
playing critical exponents at a dynamic phase transition 
which are identical to those found at the phase transition 
of the equilibrium Ising model 24[ . This surprising ob- 



servation is consistent with a symmetry argument given 
in [30(| and has been substantiated through the study of 
the time-dependent Ginzburg-Landau model in an oscil- 
lating field [31j . 

In the past very few studies have looked at the im- 
pact surfaces can have on this dynamic phase transi- 
tion, thereby focusing mostly on rather general aspects. 
For example, the effects of boundaries on magnetization 
switching in kinetic Ising models were studied in [32| ■ In 
"ID, 33| the dynamic phase transition was investigated in 



Heisenberg films with competing surface fields. 

In this Letter we present the first study of the surface 
critical properties at a dynamic phase transition. Using 
large-scale numerical simulations, we study kinetic Ising 
models with free surfaces subjected to a square-wave os- 
cillating field. Both in two and three space dimensions 
we obtain values for the surface critical exponents that 
differ markedly from the values of the equilibrium surface 
exponents, thus demonstrating that the dynamic surface 
universality class differs from that of the equilibrium sys- 
tem, even though the same universality class prevails for 
the corresponding bulk systems. In addition, we find that 
the kinetic surface phase diagram in three dimensions is 
remarkably simple and does not exhibit a special tran- 
sition point, nor a surface or extraordinary transition, 
which are all present in the equilibrium surface phase 
diagram. 

In order to study the surface critical behavior at the 
dynamic phase transition we consider square and cubic 
lattices with open boundary conditions in one direction, 
called z-direction in the following, whereas in the direc- 
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tion(s) perpendicular to the ^-direction we have periodic 
boundary conditions 34|. In this way we have in a sys- 



tem of linear extend L two surfaces, located at z = 1 and 
z = L. Every lattice site x is characterized by an Ising 
spin 5x = ±1. The Hamiltonian is given by 
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where Jf, > and J s > are ferromagnetic bulk and 
surface coupling constants. The first sum is over near- 
est neighbor sites where at most one of the sites is in 
a surface layer. The second sum, on the other hand, is 
over neighboring sites that are both in a surface layer. 
We thereby allow for different values of the coupling con- 
stants at the surface and inside the bulk. Finally, the last 
term is the magnetic field term where H (t) is a spatially 
uniform field that oscillates in time. We follow 24j and 
use a square- wave field with amplitude H . Both tem- 
perature and magnetic field strength are chosen in such 
a way that the system is in the multidroplet regime [3] : 
T = 0.8T c 2<i , H = 0.3J h for d = 2 and T = 0.8T c 3d , 
H = 0.4 J b for d = 3. Here T 2d = 2.269 ■ • • J b /k B and 
rp3d _ 4.5H5 j b /fcg are the critical temperatures of the 
two- and three-dimensional equilibrium systems. 

As the surfaces break spatial translation invariance, 
all quantities of interest depend on the distance to the 
surface. We therefore define local, i.e. layer-dependent, 
quantities. Thus we consider the layer magnetization av- 
eraged over one period of the external field (t]/ 2 is the 
half-period of the oscillating field): 
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m(t, z)dt , 



(2) 



with the time-dependent magnetization m(t, z) = 
L }_ 1 J^x Sx(t) of layer z, the sum being taken over all 
spins in that (d — l)-dimensional layer. The local order 
parameter is then given by (\Q(z)\) where (• ■ ■) indicates 
both a time average (i.e. an average over many periods) 
and a thermal average (realized in the numerical simu- 
lations through multiple independent runs with different 
random number sequences), yielding typically a total of 
500,000 periods over which the average is taken. In a 
similar way we define the layer Binder cumulant 



U(z) = l- 



{(Q(z)f) 
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(3) 



and the layer susceptibility 

X (z)=L d - 1 ({Q{zf)-(\Q(z)\) 2 ) . (4) 

In the following the surface quantities will be character- 
ized by an index s, whereas an index b will be given to 
the quantities from the middle of the sample. 



An important quantity in the study of the dynamic 
phase transition is the ratio 
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that quantifies the competition between the metastablc 
state, characterized by the metastable lifetime (t)^, and 
the oscillating magnetic field. For small values of we 
are in the dynamically ordered phase, whereas for large 
values the system is dynamically disordered. The quan- 
tity O therefore plays the same role as that played by 
temperature at an ordinary equilibrium phase transition. 
The metastable lifetime in our systems is again layer de- 
pendent, the value of (t(z)) being smaller in the sur- 
face layer than deep inside the bulk, as surface spins are 
coupled to fewer spins. As we are interested in the sur- 
face properties at the bulk phase transition, we define 
with respect to the bulk quantity (r)(,. The dynamic 
phase transition then takes place at the critical value 
6 C = 0.918 in d = 2 Q and Q c = 1.285 in d = 3. 

As shown in Figure 1 for the two-dimensional system 
with J s — J b , both the bulk and the surface order pa- 
rameters decrease rapidly when approaching the critical 
point O c from below. Concomitantly, the bulk and sur- 
face susceptibilities display peaks in the vicinity of <d c . 
Changing the system size yields system size dependencies 
(shifts of the positions of the maxima of the susceptibili- 
ties, increasing peak heights, • • •) typical for a continuous 
phase transition. 

In the infinite system, the order parameter and the sus- 
ceptibility show an algebraic behavior close to the critical 
point: 



(\Q b \) ~ (e c - ef , X b~|0c-e|- 



(0) 



Using finite-size scaling [35J, the authors of [24j found 
in two dimensions the same values for the exponents as 
those obtained for the equilibrium Ising model, namely 
/3/v = 1/8 and j/u — 7/4, where v is the critical ex- 
ponent that governs the divergence of the correlation 
length. 

Similarly, surface critical exponents are introduced to 
describe the behavior of surface quantities close to the 
bulk critical point (we use here the standard nomencla- 
ture of surface critical phenomena, see [34, 3(1 37 1): 



(|Q a |>~(e c -e)^ , xs-iee-er 11 



(7) 



Close to a bulk critical point, finite-size scaling theory 
34, 36, 37] provides us with scaling relations for our sur- 
face quantities: 



k\Qs\) = l-W v f ± {olV v ) 



(8) 

XI = h'^>~u ± {Vh~>-) (9) 

where 9 = ^ e e e< ^ , whereas F± and G± are scaling func- 
tions, where the + (— ) sign corresponds to > (<) C . 
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FIG. 1: (Color online) Bulk (a,b) and surface (c,d) quantities 
for the two-dimensional model, composed of L x L spins, with 
J a = Jb- Close to the bulk critical point C = 0.918, the local 
order parameters decrease rapidly and the local susceptibili- 
ties display pronounced maxima. Here and in the following 
error bars are smaller than the sizes of the symbols. 



Choosing = C , we therefore expect that our quanti- 
ties depend algebraically on the linear system size. This 
is shown in Figure [2] for various values of the surface cou- 
pling J s . For not too large values of J s corrections to 
scaling are negligible, so that we can determine the val- 
ues of the critical exponents from the slopes. We find 
Pl/v = 0.43(1) and m/v = 0.18(1). We immediately 
remark that these values rather well fulfill the scaling 
relation 2/3\ + 7n = d — 1 that is expected to hold for 
surface critical exponents. We also note that our values 
differ strongly from the values of the surface exponents 
in the equilibrium critical Ising model: /3f/u = 0.5 and 



7f 1 /f = [34, 36, 37]. We therefore have the interesting 
situation that while the dynamic phase transition in the 
bulk belongs to the universality of the equilibrium Ising 
model, this is not true for the corresponding surface uni- 
versality class. 

For the three-dimensional system, the bulk system un- 
dergoes again a dynamic phase transition characterized 
by the critical exponents of the three-dimensional equilib- 
rium critical Ising model, as we verified. For the surface, 
however, the situation is more complicated. For not too 
small values of the surface coupling, the situation is simi- 
lar to the two-dimensional case, see the example J s = 2Jb 
shown in the first row in Fig. [3j When reducing the 
surface undergoes at the bulk transition value C a tran- 
sition to a dynamically ordered phase. This transition is 
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FIG. 2: (Color online) Log-log plot of (a) the surface order 
parameter and (b) the surface susceptibility as a function of 
the linear system size for the two-dimensional kinetic Ising 
model at = C . The different curves correspond to different 
values of the surface coupling constant J 3 (given in units of 
Jb). The dashed lines have slopes —0.43 (a) and 0.18 (b). For 
large values of J s corrections to scaling become sizeable. 



revealed by a characteristic peak in the surface suscepti- 
bility as well as by a crossing at C of the surface Binder 
cumulant computed for different system sizes. However, 
for values of J s < 1.5J&, see Fig. [3b, the surface spins do 
not order dynamically at C , but instead are still able of 
following the external field, even though this is no longer 
the case for the bulk spins. At lower values of the sur- 
face order parameter of our finite systems deviates from 
zero, but this partial dynamical ordering is not related 
to a phase transition. This is also revealed by the pres- 
ence of a non-critical peak (or, for very small values of 
J s , by the complete absence of any peak) in the surface 
susceptibility (see Fig. [3bl) as well as by the absence of 
the crossing of the surface Binder cumulants for different 
system sizes, as shown in the inset of Fig. |3Jd. 

Based on our data, where we studied surface couplings 
from J s = to J s = 16Jb, the surface phase diagram 
of the three-dimensional kinetic Ising model shown in 
Fig. 2Jd differs remarkably from the corresponding di- 



agram of the equilibrium model [34], 1361 . 137} . see Fig. 
I3K- Not only does the surface not order dynamically at 
the bulk critical point for J s < 1.5J&, as just discussed, 
the kinetic Ising model also does not exhibit a surface 
transition, where the surface orders alone, whereas the 
bulk remains disorder. Concomitantly, the special tran- 
sition point, where both surface and bulk are critical, 
and the extraordinary transition, where the bulk orders 
in presence of an ordered surface, are also absent. In fact, 
whereas for the equilibrium system is it possible to shift 
the phase transition temperature of the two-dimensional 
surface /cbT c m = 2.269 ■ • • J s above the bulk transition 
temperature hsT^ — 4.51 15 J b by sufficiently increas- 
ing the ratio J s /Jb of the couplings, a similar mechanism 
does not exist in the kinetic Ising model. 

Finally, for J s > 1.5J& we can again measure the 
surface critical exponents through a finite-size scaling 
analysis. As shown in Fig. [5] corrections to scaling 



4 




1.6 0.8 1 1.2 1.4 




6 0.8 1 1.2 1.4 




FIG. 3: (Color online) Surface susceptibility (a) and surface 
Binder cumulant (b) for the three-dimensional kinetic Ising 
model with J a = 2J D . For small values of J a , the surface does 
not order dynamically at the bulk critical point, as shown by 
(c) the surface order parameter and (d) the surface suscepti- 
bility. All surface coupling constants are expressed in units of 
Jb- The inset in (b) shows that for J s — Jb the surface Binder 
cumulants do not cross at a common value of 0. The system 
size in (c) and (d) is L — 96. 
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FIG. 4: (Color online) Surface phase diagram of (a) the equi- 
librium three-dimensional Ising model and (b) the nonequi- 
librium three-dimensional kinetic Ising model. 



are much more important in three than in two dimen- 
sions. Based on our data, we obtain fti/v — 0.88(3) and 
ln/v = 0.29(3). These values again differ markedly from 
the known values Pl/v = 1.27 and 7^ = —0.40 of the 
corresponding surface critical exponents [38|. Most no- 
tably, whereas in the equilibrium system the surface sus- 
ceptibility displays a cusp singularity characterized by 
a negative critical exponent, in our system the surface 
susceptibility diverges with a positive critical exponent. 

As we mentioned in the introductory remarks, a sym- 
metry argument put forward in (30j states that contin- 




30 



20 



1 

• J s = 2 


(b) 


□ J =16 




1 





50 



100 



FIG. 5: (a) Surface order parameter and (b) surface suscepti- 
bility as a function of system size at the bulk critical point of 
the three-dimensional kinetic Ising model. The dashed lines 
have slopes -0.88 (a) and 0.29 (b). 



uous transitions with up-down symmetry and noncon- 
served order parameter should fall into the university 
class of the ordinary Ising model. This indeed agrees 
with our own results (as well as with previous results 

2J, |3l| ) that the bulk critical exponents at the dynamic 
phase transition are the same as that of the equilibrium 
Ising model, and this both in two and three dimensions. 
However, once surfaces are introduced, the lattice sym- 
metry is broken close to the surfaces, and one of the 
assumptions underlying the argument of [30| is no longer 
fulfilled. Indeed, our results show that the dynamic sur- 
face exponents differ from the surface exponents of the 
equilibrium model, yielding new nonequilibrium surface 
universality classes. Using field-theoretical methods sim- 
ilar to those developed for equilibrium critical surfaces 

37 , [38| , it should be possible to compute these new ex- 
ponents and to classify the possible dynamic surface uni- 
versality classes. 

It follows from our work that our understanding of the 
role played by surfaces in nonequilibrium systems, and 
more specifically at nonequilibrium phase transitions, is 
far from being complete. Surfaces break lattice symme- 
tries, and this can have many surprising and unexpected 
effects out of equilibrium, as exemplified in our study of 
surface critical behavior at a dynamic phase transition. 
Based on our results we expect that future in-depth stud- 
ies of the role of surfaces far from equilibrium will reveal 
additional new and unexpected phenomena. 

This work was supported by the US National Science 
Foundation through DMR-0904999. 



[1] T. Chou, K. Mallick, and R. K. P. Zia, Rep. Prog. Phys. 
74, 116601 (2011). 

[2] B. Schmittmann and R. K. P. Zia, Statistical Mechan- 
ics of Driven Diffusive Systems, Phase Transitions and 
Critical Phenomena , Vol. 17, eds. C. Domb and J. L. 
Lebowitz (Academic Press, New York, 1995). 

[3] H.-K. Janssen, B. Schaub, and B. Schmittmann, Z. Phys. 



■5 



B 71, 377 (1988). [23 
R. Harris and M. Grant, Phys. Rev. B 38, 9323 (1988). 
J. W. Essam, A. J. Guttmann, I. Jensen, and D. Tan- [24 
laKishani, J. Phys. A 29, 1619 (1996). 

P. Frojdh, M. Howard, and K. B. Lauritsen, Int. J. Mod. [25 
Phys. B 15, 1761 (2001). 

M. Pleimling and F. Igloi, Phys. Rev. Lett. 92, 145701 [26 

(2004) . [27 
M. Pleimling, Phys. Rev. B 70, 104401 (2004). 

M. Pleimling and F. Igloi, Phys. Rev. B 71, 094424 [28 

(2005) . 

F. Baumann and M. Pleimling, Phys. Rev. B 76, 104422 [29 
(2007). 

D. Kadau, A. Hucht, and D. E. Wolf, Phys. Rev. Lett. 

101, 137205 (2008). [30 

A. Hucht, Phys. Rev. E 80, 061138 (2009). 

F. Igloi, M. Pleimling, and L. Turban, Phys. Rev. E 83, [31 
041110 (2011). 

E. Agliari, M. Casartelli, and A. Vezzani, Eur. Phys. J. [32 
B 60, 499 (2007). 

E. Agliari, M. Casartelli, and A. Vezzani, Eur. Phys. J. 
B 84, 317 (2011). [33 
M. Pleimling, B. Schmittmann, and R. K. P. Zia, EPL 
89, 50001 (2010). [34" 
L. Li and M. Pleimling, EPL 98, 30004 (2012). [35 
M. Henkel, H. Hinrichsen, and S. Liibeck, Non- 
Equilibrium Phase Transitions: Volume 1: Absorb- 
ing Phase Transitions (Springer/Dordrecht and Cano- [36 
pus/Bristol, 2008). 

G. Odor, Universality in Nonequilibrium Lattice Systems: [37 
Theoretical Foundations (World Scientific/Singapore, 
2008). 

B. Chakrabarti and M. Acharyya, Rev. Mod. Phys. 71, [38 
847 (1999). 

M. Acharyya, Int. J. Mod. Phys. C 16, 1631 (2005). 
M. Acharrya, Phys. Rev. E 56, 2407 (1997). 



S. W. Sides, P. A. Rikvold, and M. A. Novotny, Phys. 
Rev. Lett. 81, 834 (1998). 

G. Korniss, C. J. White, P. A. Rikvold, and M. A. 
Novotny, Phys. Rev. E 63, 016120 (2000). 

H. Jang, M. J. Crimson, and C. K. Hall, Phys. Rev. B 
67, 094411 (2003). 

M. Acharyya, Int. J. Mod. Phys. C 14, 49 (2003). 

O. Canko, B. Deviren, and M. Keskin, J. Phys. Cond. 

Matt. 18, 6635 (2006). 

Q. Jiang, H.-N. Yang, and G.-C. Wang, Phys. Rev. B 52, 
14911 (1995). 

D. T. Robb, Y. H. Xu, O. Hellwig, J. McCord, A. Berger, 
M. A. Novotny, and P. A. Rikvold, Phys. Rev. B 78, 
134422 (2008). 

G. Grinstein, C. Jayaprakash, and Y. He, Phys. Rev. 
Lett. 55, 2527 (1985). 

H. Fujisaka, H. Tutu, and P. A. Rikvold, Phys. Rev. E 
63, 036109 (2001). 

H. L. Richards, M. Kolesik, P.-A. Lindg0ard, P. A. 
Rikvold, and M. A. Novotny, Phys. Rev. B 55, 11521 
(1997) 

H. Jang and M. J. Crimson, Phys. Rev. E 63, 066119 
(2001). 

M. Pleimling, J. Phys. A: Math. Gen. 37, R79 (2004). 
K. Binder, in Finite-Size Scaling and Numerical Simula- 
tions of Statistical Systems, edited by V. Privman (World 
Scientific, Singapore, 1990), p. 173. 

K. Binder, in Phase Transitions and Critical Phenomena 
Volume 8 (London/New York, Academic Press, 1983). 
H. W. Diehl, in Phase Transitions and Critical Phe- 
nomena Volume 10 (London/New York, Academic Press, 
1986). 

H. W. Diehl and M. Shpot, Nucl. Phys. B528, 595 
(1998). 



